Quasi-Equatorial Gravitational Lensing by Spinning Black Holes 

in the Strong Field Limit 
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Spherically symmetric black holes produce, by strong field lensing, two infinite series of rela- 
tivistic images, formed by light rays winding around the black hole at distances comparable to the 
gravitational radius. In this paper, we address the relevance of the black hole spin for the strong 
field lensing phenomenology, focusing on trajectories close to the equatorial plane for simplicity. In 
this approximation, we derive a two-dimensional lens equation and formulae for the position and 
the magnification of the relativistic images in the strong field limit. The most outstanding effect 
is the generation of a non trivial caustic structure. Caustics drift away from the optical axis and 
acquire finite extension. For a high enough black hole spin, depending on the source extension, we 
can practically observe only one image rather than two infinite series of relativistic images. In this 
regime, additional non equatorial images may play an important role in the phenomenology. 

PACS numbers: 95.30.Sf, 04.70.Bw, 98.62.Sb 
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I. INTRODUCTION 

Gravitational lensing, since its beginning, has been 
used to test General Relativity in the weak field approx- 
imation. Later on, several studies of null geodesies in 
strong gravitational fields have been lead in past years. 
Bardeen et al. 1] studied the appearance of a black hole 
in front of a uniform source; Viergutz 0] made a semi- 
analytical investigation of null geodesies in Kerr geome- 
try; Nemiroff |3| studied the visual distortions around a 
neutron star and around a black hole; Falcke, Melia & 
Agol Q| considered the luminosity of the accretion flow 
as a source. 

A recent paper by Virbhadra & Ellis [5j has risen a new 
interest about gravitational lensing as a probe for strong 
gravitational fields generated by collapsed objects, pro- 
viding a new important test for the full general relativity. 
They have shown that a source behind a Schwarzschild 
black hole would generate an infinite series of images 
on both sides of the lens. These relativistic images are 
formed by light rays passing close to the event horizon 
and winding several times around the black hole before 
emerging towards the observer. By an alternative formu- 
lation, Frittelli, Kling & Newman [(| attained an exact 
lens equation, giving integral expressions for its solutions, 
and compared their results to those by Virbhadra & El- 
lis. The phenomenology of collapsed objects with naked 
singularities, analyzed by Virbhadra & Ellis in another 
work |2j, is radically different. This difference provides a 
possible way to test the correctness of the cosmic censor- 
ship conjecture. 

A new, simple and reliable method to investigate the 
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subject was proposed by Bozza et al. in Ref. They 
revisited the Schwarzschild black hole lensing defining a 
strong field limit for the deflection angle, which retained 
the first two leading order terms. By this approxima- 
tion, a fully analytical treatment was developed and sim- 
ple formulae for the position and the magnification of 
the images were derived. The same method was applied 
by Eiroa, Romero & Torres ^ to a Reissner-Nordstrom 
black hole, confirming the appearance of a similar pat- 
tern of images. Later on, the formulae given in Ref. || 
were used by Petters 01 to calculate relativistic effects 
on microlensing events. Finally, in a previous work, we 
have developed a generalization of the strong field limit to 
an arbitrary spherically symmetric spacetime com- 
paring the image patterns for several interesting metrics. 
We have shown that different collapsed objects are dis- 
tinguishable by a careful examination of the separation 
between the first two relativistic images and their lumi- 
nosity ratio. The sensitivities required for such measures 
are out of the actual VLBf projects .12], but might be 
reached in a not so far future. 

Insofar, only spherically symmetric black holes have 
been adequately investigated. Yet, in general, a black 
hole would be characterized by a non-zero intrinsic angu- 
lar momentum which breaks spherical symmetry, leaving 
only a rotational symmetry around one axis. Rotation 
heavily affects the gravitational field around the collapsed 
object. It is thus natural to expect relevant modifications 
in the phenomenology of strong field gravitational lens- 
ing. 

A further motivation for such a study comes from the 
fact that previous works on the subject have selected the 
supermassive black hole hosted by the radio source Sagit- 
tarius A* 0] as the best candidate for strong field gravi- 
tational lensing §,[n|. Its mass has been estimated to be 
M = 2.6 x 10 6 M Q but our knowledge on a possible intrin- 
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sic angular momentum of this object is still very poor. 
However, the analysis of the variability of the spectrum 
in the mm-submm reg ion suggests the possibility of a 
non-negligible spin |14| . In particular, the value 



0.044 



(1) 



has been proposed (with respect to Ref. [Ij] , we normal- 
ize distances to the Schwarzschild radius rather than to 
the black hole mass, hence the factor 1/2 in this value), 
but high uncertainties in the assumptions behind the cal- 
culations may push the spin towards even higher values. 

The purpose of this paper is to investigate the rele- 
vance of the black hole spin in strong field lensing phe- 
nomenology. We formulate the strong field limit for Kerr 
black hole lensing, considering trajectories close to the 
equatorial plane. With this limitation we cannot give a 
complete description of the whole phenomenology, which 
we shall delay to future works. However, as we shall see, 
quasi-equatorial motion yields very sharp indications use- 
ful to understand also the general case. In particular, we 
will show that the presence of a sufficiently high black 
hole spin drastically changes the expected pattern of ob- 
servable images. 

This paper is structured as follows. In Sect. 2 we re- 
call some general properties of Kerr geodesies. In Sect. 
3 we carry out the strong field limit expansion of the 
deflection angle on the equatorial plane. In Sect. 4 we 
move off from the equatorial plane and consider trajec- 
tories at small declinations. In Sect. 5 we write the lens 
equation on the equatorial plane. In Sect. 6 we write a 
polar lens equation dealing with displacements normal to 
the equatorial plane. In Sect. 7 we find the positions of 
the caustic points and the magnification of all images in 
the equatorial plane. In Sect. 8 we describe the critical 
curves and caustic structure. Sect. 9 discusses the effects 
of the black hole spin on the gravitational lensing phe- 
nomenology, with special reference to the black hole at 
the center of our Galaxy. Sect. 10 contains the summary. 



II. GEODESICS IN KERR SPACETIME 

In Boyer-Lindquist coordinates 0] x^ = (t,x,i}, 
the Kerr metric reads 



ds 2 



A - a 2 sin i? , 9 p , . > , > 
- dt 2 - '--dx 2 - p 2 dd 2 



2 i „2\* r ,2 A 



(X + a 



a 2 A sin 2 # 



■ sin 2 $d(j> 2 



A = x 2 - x + a 2 
p 2 = x 2 + a 2 cos 2 ■& 



2ax sin 



dtd(f) (2) 

(3) 
(4) 



where a is the specific angular momentum of the black 
hole. All distances are measured in Schwarzschild radii 



The Kerr space is characterized by a spherical event 
horizon at xh 



l + \/l-4« 2 



for \a\ < 0.5. Beyond this 
critical value of the spin there is no event horizon and 
causality violations are present in the whole spacetimc, 
with the appearance of a naked singularity [Hi [l3| • We 
shall restrict to subcritical angular momenta. The ellip- 
soid A — a 2 sin 2 1? = is the static limit bounding the 
region where every static worldine (x = $ — 4> — 0) 
becomes spacelike. The region between the static limit 
and the horizon is called ergosphere: here everything is 
bound to rotate around the black hole. 

The geodesies equations can be derived taking 



£ — Qiii/X^ 



(5) 



as the Lagrangian, where the dot indicates the derivative 
with respect to some affine parameter. 

Finding four integrals of motion, we can transform 
these equations into a set of four first order equations 
which are equivalent to the original ones. Two costants 
of motions are the energy and the angular momentum of 
the particle, given by 



(0) 
(7) 



By a suitable choice of the affine parameter, we can set 

E=l. (8) 

From these equations, we find an expression for t and 
6 in terms of x, 1) and J 



333 + gozJ 

333500 - .9o3 
ffOO J + ff03 
503 ~~ ff33fftt 



(9) 
(10) 



C is another constant of motion, which vanishes for 
null geodesies and can be used to write x 



/-.9oo< 2 - 922-d 2 - .g 33 2 - 2g 3^ ,,,, 
x = ±\ . (11 

V 511 

Finally, $ can be obtained in terms of a fourth integral 
of motion, separating the Hamilton- Jacobi equation [16j : 



i) = ±\\/Q + a 2 cos 2 d-J 2 cot 2 1?. 
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1) 



(12) 



Eqs. lffi]l-l|12[l represent the sought set of first order 
differential equations, suitable for a detailed study. 

The integrals of motion J and Q can be expressed in 
terms of the geometric parameters of the incoming light 
ray trajectory. In general, we can identify a light ray 
coming from infinity by three parameters (Fig. refer- 
ring to the straight line which the photon would follow if 
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FIG. 

ray. 



1: The three parameters identifying an incoming light 



there were no gravitational field. The projection of this 
line on the equatorial plane has a distance u from the 
origin, which we shall call the projected impact param- 
eter. At this minimum projected distance, the light ray 
has some height h on the equatorial plane. Finally, the 
inclination -ipo is the angle that the light ray forms with 
the equatorial plane. When we switch on the gravita- 
tional field, the light ray is obviously deviated from this 
ideal straight line, but these three parameters can still 
be used to label any light rays coming from infinity. 

Taking the asymptotic limit of the equations of motion, 
we can write J and Q in terms of the initial conditions 
characterizing the light ray 



III. 



J = u cos ipo 

Q = h 2 cos 2 tpo + (u 2 - a 2 ) sin 2 ip . 



DEFLECTION ANGLE IN THE 
EQUATORIAL PLANE 



(13) 
(14) 



In this section, we consider light rays strictly lying on 
the equatorial plane d — ? by setting h = i/jq = 0. The 
reduced metric has the form 

ds 2 = A(x)dt 2 ~ B(x)dx 2 - C(x)d(/) 2 + D(x)dtdcf>. (15) 

with 

A(x) = 1-1 
B(x) = - 



-73- 



(16) 



C(x) = x 2 + a 2 + ^ 
D(x)=2%. 

What we say in this section is immediately extendable 
to any axially symmetric spacetime if we replace Eq. 1161) 
by any other expression. 



As ipo = 0, by Eq. the angular momentum J coin- 
cides with the impact parameter u. The impact param- 
eter is also related to the minimum distance xq reached 
by the photon. In general, a light ray coming from in- 
finity approaches the black hole, reaches this minimum 
distance £0 and then leaves again towards infinity. Eval- 
uating the Lagrangian at x = xo, we find an implicit re- 
lation between J = u and the closest approach distance 
x 



J = u 



-Do + y/AApCp + Dl 
2A Q 



a + x Q ^/a 2 + x (xo - 1) 



x 



1 



(17) 



where all the metric functions with the subscript are 
evaluated at x — xq. The impact parameter u is then 
univocally determined by xq and vice versa. Choos- 
ing the positive sign before the square root, we describe 
only light rays winding counterclockwise when seen from 
above. For a > the black hole also rotates counter- 
clockwise, while for a < the black hole and the photon 
rotate in opposite senses. 

Dividing <j) by x, we find the azimuthal shift as a func- 
tion of the distance 



d(j) x 

— = Pi{x,Xo)P2(X,Xq) 

dx 



Pi(x,x ) 



P2{x,x Q ) 



B(2A AJ + A D) 

Vca^VIacTd^ 

1 



A -A^ + ^(AD Q ~A D) 



(18) 
(19) 
(20) 



Integrating this expression from xq to infinity we find 
half the deflection angle as a function of the closest ap- 
proach. Given the symmetry between approach and de- 
parture, we can write the whole deflection angle as 



a(x ) = <j>f(x ) - 7r 



4>f{xo) 



2 / --dx. 

dx 



(21) 
(22) 



4>f{xo) is the total azimuthal shift. It evaluates to 7r 
for a straight line and becomes larger as the light ray 
is bent by the gravitational field. The expression for a 
spherically symmetric metric, given in Ref. can be 
recovered setting D — Dq — 0. 

The deflection angle grows as xo decreases. It diverges 
when xq reaches a minimum value x m which represents 
the radius of the photon sphere. If a photon falls inside 
this sphere, it is destined to be absorbed by the black 
hole. Of course, we will have different photon spheres 
for photons winding in the same sense of the rotation of 
the black hole (hereafter direct photons) and for photons 
winding in the opposite sense (retrograde photons). In 
general, we expect the latter to be absorbed more easily. 
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Their photon sphere will thus be larger than that of ret- 
rograde photons, which can escape more easily. As we 
shall see later, this is what happens. 

Following the philosophy of the strong field limit, we 
look for an expansion of the deflection angle of the form 

a(9) = -a log (0E°L - l ) +b + 0{u-u m ) (23) 

V Um J 

where the coefficients u m , a and b depend on the metric 
functions evaluated at x m . Dol is the distance between 
the lens and the observer, so that the angular separation 
of the image from the lens (also referred as impact angle) 

fefl = _!!_. 

L>OL 

All the steps to be taken towards this final expression 
are very similar to those for spherically symmetric black 
holes, with few adjustments. We shall sketch them very 
briefly, referring to Ref. ^lj f° r details. 

We define the variables 



y 



A(x) 

y-yo 

l-l/o 



(24) 
(25) 



where y$ = Aq. The integral (|22|1 in the deflection angle 
becomes 



<f>f(xo) = J R(z,x )f(z,x )dz (26) 
o 

R(z,x ) = 2±—^P 1 (x 1 x ) (27) 

f(z,x )=P 2 {x,xo) (28) 

where x = A^ 1 [(1 — yo) z + yo\- 

The function R(z, Xq) is regular for all values of z and 
xq, while f(z,xo) diverges for z — > 0. To find out the 
order of divergence of the integrand, we expand the ar- 
gument of the square root in f(z, Xq) to the second order 




a 



FIG. 2: The radius of the photon sphere versus the black hole 
angular momentum. 



The real solution external to the horizon of this equation 
defines the radius of the photon sphere x m , plotted in 
Fig. As expected, for positive a, (direct) photons are 
allowed to get closer to the black hole, entering even the 
ergosphere (x m falls below 1) at high values of a. It is 
possible to calculate exactly from Eq. a = at what 
angular momentum this happens. The critical value is 



1 



2V2 



= 0.354. 



(32) 



Notice that x m — > 1/2 as a — ► 1/2, i.e. the photon 
sphere coincides with the horizon in the limit of extremal 
Kerr black hole. For negative angular momenta, (retro- 
grade) light rays must keep farther from the center. 

The procedure to find the strong field limit coefficients 
is from now on identical to that described in Ref. [Tlj . 
with R(z,xq) f(z,Xo) and f(z,Xo) given by Eqs. J57J), 
(|2*8|) and , respectively. We shall not repeat the whole 
technique here but just specify the results of Ref. 01 f° r 
our metric. 

The strong field limit coefficients of the expansion ll'MI) 



f(z,x ) ~ fo{z,x ) 



1 



+ (iz 2 



(29) 



When a is non zero, the leading order of the divergence 
in fo is z -1 / 2 , which can be integrated to give a finite 



result. When a vanishes, the divergence is 



which 



makes the integral diverge. Then the outermost solution 
of the Eq. a = defines the radius of the photon sphere 
x m (see also 0). 

In the case of the Kerr metric, we have 



a = Xq 



where 



:ro(3 — 5xo + 2xq) — 2a 2 + 2a^J a 2 + xq{xq — 1) 
•[(x -l)(x3 + a 2 (x + l)] _1 . (30) 
Eq. a = is equivalent to the third degree equation 



-D rn + y/AA m C m + £> 2 Tl 



2A„ 



_ R(0,x m ) 

a — 



-7r + 6 D + bfl + alog— ^ 



bo = 2a log 



2(1 -y m ) 
A' m x m 



(33) 
(34) 
(35) 

(36) 



b R = I [R{z,x m )f(z,x m ) - R(0,x m )f o (z,x m )] d$7) 



and c is defined by the expansion 



8a z 



x (3 - 2x ) 2 = 0. 



(31) 



M-«ra = c(l - X m y 



(38) 
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IV. PRECESSION OF THE ORBITS AT SMALL 
DECLINATIONS 
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FIG. 3: Coefficients of the strong field limit versus black hole 
angular momentum. 



All the functions with the subscript m are evaluated at 

Xq — X m . 

Fig. |31 shows the strong field limit coefficients as func- 
tions of a. The minimum impact parameter u m decreases 
with a in a way similar to x m . a grows and b decreases, 
both diverging with opposite signs at a = 1/2. The diver- 
gence of the coefficients of the expansion warns that the 
strong field limit deflection angle i12-'il) no longer repre- 
sents a reliable description in the regime of high a, since 
x m becomes a higher order singularity in the function 
f(z,xo) and a different expansion should be performed. 



The study of the deflection angle for photons lying 
on the equatorial plane is sufficient to write a one- 
dimensional lens equation. However, to address the caus- 
tic structure and calculate the magnification of the im- 
ages, we need a two-dimensional lens equation. For this 
reason, in this section we analyze trajectories close to 
the equatorial plane. They are described by one further 
coordinate: the polar angle i?, or equivalently, the dec- 
lination if) = f — 1?- The problem becomes too involved 
to be solved in general but we shall give a complete de- 
scription for quasi-equatorial motion, preserving also the 
simplicity and immediacy of the strong field limit scheme. 

In order to remain at small declinations, we restrict 
to light rays characterized by a small inward inclination 
?po and small height h compared to the projected impact 
parameter u, with ipQ ~ — . Retaining the first relevant 
terms, from Eqs. {I3J|-|(TI$ we get 



J '. 

Q 



—2 / 2 



U — \J u 2 



(39) 
(40) 
(41) 



We require the declination ip to stay small (of the order 
of tpo) during the motion. Dividing Eq. I|12|) by Eq. I|1(J|) . 
we get a simple evolution equation for ip as a function of 
the azimuth d> 



with 



w(0) 



h2 ,i 

=2 + ^0 



_ a 2 + x(<f>)(x(4>) - 1) 



(42) 

(43) 
(44) 



In the Schwarzschild case (a 
is immediately solved to 



a + u{x((f>) - l))x(4>)' 

0), uj -> 1 and Eq. |g2J) 



ll){4>) =1pCOS((f) + <f>o). 



(45) 



After each loop around the black hole, the declination 
returns to the initial value. This means that there is no 
precession of the orbital plane, as expected for a spheri- 
cally symmetric black hole. 

For non-vanishing angular momenta, u> is no longer a 
constant, since x depends on (f>, and the solution of Eq. 
(|42|1 is generalized to 



with 



ip((p) = V'cos ((f) + <f> ) 



;(0W- 



(46) 



(47) 
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Since we are interested to gravitational lensing, the 
photon comes from infinity and returns to infinity. There- 
fore, we are interested to the quantity 



<Pf 



;(/)#'■ 



(48) 



where 0/ is the total azimuthal shift experienced by the 
photon in its whole trajectory, given by Eq. 1)22)1, 
This integral can be rewritten as 



= 2 



uj{x) — dx = 
dx 



Ruj(z,x )f{z,x )dz (49) 



.to 



where 



Ru(z, xq) = u(x)R(z, xo) 



(50) 



with R(z,xq) and f(z,xo) given by Eqs. 1)27(1 and 1)28( 1. 
respectively. At this point, the integral can be solved 
by the same technique used for the integral ((26)1 with a 
very similar result. This is possible since u>{x) adds no 
singularities and can be englobed in the regular function 
R. 

The final result is 



4>f = -a log I 1 



Ruj (0; %r. 



= 1 



-7T + b 



D i- Or- 



aloj 



where 



&d = 2a log 



2(1 - y m ) 

A m x m 



(51) 
(52) 
(53) 

(54) 



b R = J [R u (z,x m )f(z,x m ) - R u (0,x m )f (z,x m )]pi5) 
o 

Notice that the coefficient of the logarithmic term 
a turns ro be exactly equal to 1 for all values of the 
black hole spin. The coefficient b is just b + tt in the 
Schwarzwschild limit a — 0. This recovers the equiva- 
lence between the phase in the polar motion </>j and the 
total azimuthal shift 0/ in this limit. The full behaviour 
of b is plotted in Fig. This coefficient diverges as well 
in the extremal black hole limit a — > 1/2. 

As a general remark, we can say that for positive angu- 
lar momenta lo is always less than one, so that (j>f < <frf. 
As a consequence, the orbital plane suffers a counter- 
clockwise precession, i.e. after each loop it is necessary 
an additional A</> to reach the same declination ip. On the 
contrary, for negative angular momenta, lo > 1. In this 
case, the precession is clockwise, i.e. the photon reaches 
the same declination before completing a loop. 



b 
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FIG. 4: The coefficient b as a function of a. 



The integration constant 4>q in Eq. 1)46)1 is fixed by the 
initial conditions. In particular, we have to impose that 
at (j) — the declination is just minus the inclination of 
the incoming photon trajectory, that we have indicated 
by if>Q. The result is that 



a 



-Sign[/i] arccos 



(56) 



The declination of the outward photon is thus 

V>/ = = ?cos fa + O ) ■ (57) 

In alternative, using the expression of 4>o, we can write 



h . - 

ijjr = —ipo cos Of — — sin ( 



7- 



(58) 



The phase <fi f , calculated through Eq. 1)51)1, has a cen- 
tral importance in the discussion of Sects. IVII IVIII 



V. LENSING IN THE EQUATORIAL PLANE 

In the previous sections we have expressed the deflec- 
tion angle a as a function of the impact angle 8 and the 
outward declination -0/ as a function of the incoming in- 
clination ipo- We are thus ready to write a lens equation 
for the Kerr black hole. In this section we shall write the 
equatorial lens equation, while the next section will deal 
with the polar lens equation. 

So let us start from the ideal case when observer and 
source both lie on the equatorial plane of the Kerr black 
hole and the whole trajectory of the photon is confined on 
the same plane. In previous works, the strong field limit 
has been developed assuming an almost perfect align- 
ment of source, lens and observer. This because, for 
spherically symmetric metrics, the better the alignment 
the higher is the magnification. As we shall see in Sec. 
rvTTI this is no longer the case for Kerr black holes. There- 
fore we shall write the equatorial lens equation in a more 
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FIG. 5: The lensing geometry projected on the equatorial 
plane in the case of the first relativistic image. 9 is the impact 
angle as seen by the observer, is the angular position of the 
source as seen by the observer, 7 is the angular position of 
the source as seen by the lens, 9 is the impact angle from the 
source. 



general way, allowing for a generic geometric disposition 
of lens, source and observer (see Fig. EJ. 

The optical axis is the line joining the observer and 
the lens. Setting the origin on the black hole, the angle 
between the direction of the source and the optical axis 
will be indicated by 7. 7 ~ is the case of almost perfect 
alignment discussed in Rcfs. II II 113 EH Fr om the 
lensing geometry, illustrated in Fig. we can write the 
relation 



where 



6 + 6 mod 2tt, 



Dql, 



Dls D 



LS 



(59) 



(60) 



is the impact angle from the source and Dls is the dis- 
tance between the lens and the source. 
The equatorial lens equation is then 



D 



7 : 



OL 



D 



LS , 



D LS 



— a(6) mod 2tt 



(61) 



In this lens equation 7 can assume any value in the 
trigonometric interval [— tt,tt]. The source may even be 
on the same side of the observer when 7 = n. The rela- 
tion between 7 and (3 (the angular position of the source 
as seen by the observer) is 



. R D LS . 
Bin p = — — sm 7, 
Dos 



(62) 



but here in general we cannot substitute the sines by 
their arguments. -Dos is the distance between source 
and observer which does not coincide with the distance 
Dol + Dls actually covered by lensed photons. 

To solve the lens equation, since 6 = -jj^ < 1, in a 
first step we solve the Eq. 7 = —a(6) mod 27r. Using the 



expression for the deflection angle derived in the strong 
field limit (J23J, we find 



Dol 

b+-,-2 



(l + e„) 



(63) 
(64) 



where n = 1,2,... indicates the number of loops done by 
the photon around the black hole. This solution is then 
corrected expanding a(6) around 



(6-6^ + 0(6-6^) 



-J- 



(65) 



Substituting in l|61l) and neglecting higher order terms, 
we find 



,j„ f 1 _ u m e n (D OL + Dls) ^ ^ ^ 



clDolDls 



where the correction is much smaller than Q®. 

Images are formed on both sides of the lens. As all 
strong field limit coefficients depend on a we have to be 
careful and choose the correct sign for the angular mo- 
mentum. Conventionally we call north the direction of 
the black hole spin. Then photons winding counterclock- 
wise are direct and are described by a positive a. They 
form images on the eastern side of the black hole. Images 
formed by retrograde rays appear on the western side and 
are described taking a negative a and reversing the sign 
of 7. 



VI. LENSING AT SMALL DECLINATIONS 

The lens equation l|61|l describes trajectories lying on 
the equatorial plane and can be employed to calculate 
the positions of the relativistic images. However, to in- 
vestigate the problem on a deeper level we are forced 
to study what happens at least for small displacements 
from the equatorial plane. In this section we shall assist 
Eq. lffTT|l by its polar counterpart, which is necessary to 
understand the caustic structure and compute the mag- 
nification of the images. 

Thanks to the small declination hypothesis, at the low- 
est order we can neglect any backreaction on the equa- 
torial lens equation. In all our discussion we shall speak 
(using time reversal) as the photon were emitted by the 
observer and absorbed by the source. 

Consider a source whose height on the equatorial plane 
is hs- The height of the observer will be indicated by ho- 
We shall assume that the following hierarchy of distances 
holds (see Fig. 



u < (ho,h s ) < (Dol, D LS )- 



(67) 
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FIG. 6: The lensing geometry projected on the vertical plane. 
In this picture we have assumed 7 = for simplicity. 



Recalling the meaning of the parameters tpo and h used 
insofar to identify the incoming light ray, we can write 
down the simple geometric relation 



Dol^o 



(68) 



A similar relation holds between the outgoing photon 
parameters hf,ipf and the source position 



h s = h f + D LS i!f. 



(69) 



Given the positions of source and observer, our pur- 
pose is to determine ipo, the inclination under which the 
observer emits (sees) the light ray. 

By symmetry between the outgoing and the incoming 
parameters, Eq. 143(1 for ip can be written substituting 
ipo and h by ipf and hf 




^7 



(70) 



In this way, we can express hf in terms of ipf and then, 
by Eq. (|57|l , in terms of <pf and <fio 

hf = uipsin (<j)f + (f>o) . (71) 
Recalling Eq. (JSHJl, we also get 

hf = — uipo sin <pf — h cos 4>f . (72) 

Substituting in Eq. l(6"9"|) together with Eq. (JSSJ), we 



get 



hs = -ipouS - hC - D LS tpoC ■ 



D LS =S, 
u 



where 



(73) 

(74) 
(75) 

Finally, substituting h from Eq. i|68|) and discarding 
higher order terms, we obtain the lens equation in the 
polar direction 



S- 
C 



sm q 

■ COS( 



7- 



hs 



D 



LS 



s-c 



{D OL + D LS )C - E2I^Ls 



(76) 



In this equation ipa is directly related to the heights of 
the observer and the source. The solution is 



hs + h C n - ho^-Sn 



~(D OL + D LS )C n 



DqlDls 



(77) 



where S n and C n are S and C calculated for pf = 4>f, n - 
The phase (j)f n of the n-th image is the only quantity 
that needs to be calculated preliminarly. However, once 
the equatorial lens equation (|61|l is solved, we know the 
impact angle 9 n of the n-th image and then we can cal- 
culate (j)f by Eq. (|BT|l . 

As a consistency check we can see what we obtain in 
the Schwarzschild case when the photon completes just 
one loop around the black hole, exiting on the opposite 
side. In this case a = and ^ n = (2n + l)ir. We get 



hs - he 



/ = (2n+lK DqL+Dls' 



(78) 



which is the correct result for photons passing very close 
to the black hole, looping around it. 

The consistency of our approximation requires that 
j/jq < 1 and h <C u. From Eq. (|FJ5|l the height is 



h n = 



hsDpL — hoD LS C n 
(D OL + D LS )C n + D "'- D ^ S, 



(79) 



For a generic both constraints are automatically 
satisfied, since the second term in the denominators 
dominates and we have that tpo ~ ho/DoL and h ~ 
uho I ' Dol- However, in the neighborhood of 4>j = kir the 
denominators of the two expressions can vanish, making 
diverge both quantities. The equation 

K{rf) = u{D OL + D LS )C - DqlDlsS = (80) 

defines the positions of the caustic points. In the next 
section we will discuss this equation in connection with 
the magnification of the images formed by sources close 
to the caustic points (which we call enhanced images for 
simplicity) . 

Surprisingly, thanks to the dragging phenomenon, the 
quasi-equatorial hypothesis is nearly always satisfied, ex- 
cept for enhanced images. In this situation the quasi- 
equatorial motion hypothesis is satisfied only for partic- 
ular geometric configurations which keep tpo and h under 
control. 



VII. MAGNIFICATION AND CAUSTIC POINTS 

The magnification is classically defined as the ratio of 
the angular area element of the image and the corre- 
sponding angular area element of the source that the ob- 
server would see if there were no lens. The angular area 
element of the image is 



d 2 A, 



dOdipo . 



(81) 
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The distance covered by the photons is Dql + Dls 
and then the corresponding angular area element of the 



source is 



<PA f 



D LS d-ydh s 

{Dql + D LS y 



(82) 



In fact the source element in the vertical direction is 
dhs/(DoL+DLs)- In the horizontal direction, the source 
element is span by when seen from the lens which 
corresponds to an angle D L sdj/(D L + D L s) seen from 
the observer. If we want to compare the luminosity of 
a lensed image with the luminosity of the direct image 
(namely the source observed directly along Dos with- 
out lensing), the magnification is to be multiplied by the 
factor (°o,+D,s)\ 
D os 

Our lens application has the form 



7 = 7(0) 

hs = h s (d,ip ), 



(83) 
(84) 



where the dependence on 9 in the polar lens application 
is through <f> f and we have neglected the backreaction of 
ipo on 7. The ratio between d"fdhs and dOdipo is given 
by the modulus of the Jacobian determinant of the lens 
application 



\J\ = 



d"f dht 



89 dipo 



The magnification is then given by 

d 2 Aj _ (D OL + D LS f 1 



^ = 



d 2 Ac 



D LS 



\J\ 



(85) 



(86) 



By the equatorial lens equation i|tjl|) . retaining the 
dominant terms, we have 



d-j 



aDoL 



with 



(87) 



(88) 



In the following, it is convenient to encode the number 
of loops done by the photon within 7, in order to write 
more compact formulae for all the rclativistic images. So 
7 can assume any negative real value; 7 mod 2tt repre- 
sents the angular position of the source and n — [^-^1 
is the number of loops done by the photon. In concrete, 
two values of 7 differing by a multiple of 2tt represent 
the same source position with respect to the lens, but 
reached by photons performing a different number of 
loops around the lens. For example, 7 = is a source 
aligned behind the lens reached by a photon making no 
loop (weak field lensing); 7 = — 2ir is the same source be- 
hind the lens but reached by a photon making one loop; 
7 = — 47r is the same source for a photon making two 
loops, and so on. 



By the polar lens equation (|76(l . we have 

Olpo u 

Assembling everything together, we get 



(89) 



M = 



{Dql + D LS f 



Dm.D 



oluls 



a\u{D L+D LS )C 



DolD ls S\ ' 
(90) 

For a generic (j>j, \i = O (tj^J 1 but for the en- 
hanced images, /i may even diverge (formally for point- 
like sources) when the denominator of Eq. I|9U|I vanishes. 
The 7's where this happens are called caustic points. At 



the lowest order in 



Doi 



Eq. iJHSJ reduces to 



kir. 



(91) 



Combining the formula l|51|l for the phase 4>j with the 
formula for the deflection angle (|23|l and using the equa- 
torial lens equation at the lowest order 7 = —a(9), this 
equation becomes 



7 + 6 



kir. 



(92) 



The solutions of this equation determine the angular 
positions 7/c of the caustic points 



7fe = — b + a(b — kir). 



(93) 



For each k, we have one caustic point for direct pho- 
tons and one caustic point for retrograde photons, k = 1 
would describe the weak field caustic points, formed when 
the azimuthal shift is about n. To be coherent with 
our strong field limit approximation, we shall restrict our 
analysis to k > 2. 

Expanding the denominator of Eq. I|9U|) around the 
caustic points, we have 



- K'{j k ){j - 7fc(«)) 



DqlDls 



(7- 



7fc0))- 

u 

(94) 

To understand the nature of these caustic points, no- 
tice that in the Schwarzschild limit 7^ — > — (fc — l)7r 
and all the odd caustic points are aligned on the opti- 
cal axis behind the lens on consecutive Ricmann folds 
while the even ones are aligned before the lens. If the 
source is aligned behind the lens, the n-th image is given 
by photons doing n loops around the black hole. Set- 
ting 7 ~ — 2mr, the closest caustic point is 72n+i- Then 
we can recover the Schwarzschild magnification for the 
images created by a source behind the lens @ 



,Sch 



OL 



Dls? 



UU m G r , 



q I DqlDls 



DolDls 
Dps u 2 n e n (l + e n ) 
D 2 ol Dls 



(2?w + 7) 



(95) 
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FIG. 7: The angular positions of the first five relativistic 
caustic points: k — 2,3, 4, 5, 6 from below to above. When 
7 = 2m7r the source is behind the lens, when 7 = (2m + 1)tt 
the source is before the lens. 



In Fig. 0we plot the positions of the first five relativis- 
tic caustic points as functions of the black hole angular 
momentum. The first relativistic caustic point 72 is ob- 
tained when the photon turns around the black hole and 
comes back towards the observer. 72 is thus close to — it 
but is anticipated for negative a and delayed for positive 
a. 73 is behind the lens but, at large angular momenta, 
can move very far from the initial position. At high val- 
ues of the spin, the caustic points drift so much that they 
can even change their Ricmann fold. 

We can specify the magnification formula for the en- 
hanced images using Eq. I|94|l 



enh 



_ (Dql + Dls) 2 M fc (a) 



D 2 OL D LS 



I T — Tfc I 

Mfe(o) =u(lk(a))u m {a)e lh ( a ). 



(96) 
(97) 



The quantity ~p k regulates the magnification close to 
caustic points. The dependence on 7 has been extracted 
and has the typical I7 — 7/ c |~ 1 behaviour. The depen- 
dence on the astronomical distances Dol, Dls, Dos is 
negligible in ~p k at the lowest order in jj^-- So we can 
use ~p k as a measure of the magnifying power of the black 
hole for different enhanced images and different angular 
momenta. 

In Fig. Owe plot the magnifying power ~p 3 of the caus- 
tic point 73, which for a — generates the first relativis- 
tic image of a source behind the lens. The magnification 
grows for negative angular momenta while decreases for 
positive a almost linearly. The divergence in the magni- 
fication when a approaches its extremal value 1/2 should 
not be taken seriously, as the standard strong field limit 
approximation breaks down as explained in Sect. IIIII 

The shape of ~p k remains more or less the same for 
every k but, since e 7fc = e b ~ kr! , we have that 



A*fc 



0.043. 



FIG. 8: The magnifying power at the caustic point 73. 



The magnification of enhanced images falls quite rapidly 
as we let the photons make more and more loops. 



VIII. CRITICAL CURVES AND CAUSTIC 
STRUCTURE 

It is well known that the Jacobian of the Schwarzschild 
lens has an infinite series of Einstein rings 3, 5]. The first 
one is the classical weak field Einstein ring whose angular 
radius is 



IE = 



2D 



LS 



DqlDqs 



(99) 



The corresponding caustic is the point at 71 = 0. 

At small impact parameters we enter the strong field 
limit of the Schwarzschild lensing and the light rays wind 
around the black hole. The second Einstein ring is cre- 
ated by photons coming back towards the observer. The 
caustic is at 72 = — 7r. 

Decreasing u further, the light ray completes a loop 
and we have the third Einstein ring, whose caustic point 
is 73 = — 27r and is superposed on the first caustic point 
(on the second Riemann fold). 

Summing up, the Schwarzschild lens has a large weak 
field Einstein ring and an infinite series of concentric rela- 
tivistic Einstein rings, very close to the minimum impact 
angle #00. In the region bounded by the n-th ring and 
the n + 1-th one, the sign of the Jacobian is (— 1)" +1 . 

What changes when we turn on the spin of the black 
hole? As regards the first Einstein ring of the weak 
field limit, it is distorted and shifted. As a consequence, 
the caustic point turns into a finite extension diamond 
shaped caustic |l9j. 

For the critical curves in the strong field limit, we can 
calculate their intersections with the equatorial plane, 
which are 



qO, cr 

1 k 



1 



AD, 



OL 



D 



LS, 



gDolDls 



9 0,r 



Dql 



(l + e 7 J. 



(100) 
(101) 
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FIG. 9: The first six caustics of the Kerr lens for a = 0.1, 
marked by the thick lines between 7fc(— H) and 7fc(|a|). 



They are closer to the optical axis on the positive a side 
(for left-winding photons) and farther on the negative a 
side (i.e. for right- winding photons). Therefore critical 
curves are distorted and shifted towards the negative a 
side, that is the western side, if north is the direction of 
the spin. 

The caustics are no longer points but acquire a non 
vanishing extension. 7fc(— |a|) and 7fc(H) represent the 
intersections of the k-th caustic with the equatorial plane. 
As |7fe(— |a|)| < k-K and |7fc(|a|)| > k-K the caustic is 
shifted towards the western side. To visualize this sit- 
uation, in Fig. OH we have plotted the projections of the 
caustics on the equatorial plane, as seen from the north 
direction. The [71 (— \a\), 71(H)] caustic is the weak field 
one, which stays almost aligned on the optical axis, while 
the relativistic caustics drift in the clockwise direction. 
As k grows, the caustics become larger and farther from 
their initial position on the optical axis. Examining Fig. 
we notice that at high angular momenta the caustics 
may become very large, covering even several Riemann 
folds! 

At the lowest order in ip and neglecting any backre- 
action on the equatorial lens equation it is not possible 
to give a rigorous classification of the type of catastro- 
phes we encounter on the equatorial plane. However, 
the fact that the first caustic in the weak field assumes 
the typical diamond shape of quadrupole lenses suggests 
a similar picture for strong field caustics. If this is the 
case, then the caustic points 7& on the equatorial plane 
are cusps. This is consistent with the fact that if we let 7 
decrease below some "fk the corresponding image changes 
parity. This can happen only when at the critical point 
two images are formed with the same parity of the orig- 
inal image [20I l2l| . These images rapidly move in the 
vertical direction and are missed in our quasi-equatorial 
approximation. 



IX. PHENOMENOLOGICAL IMPLICATIONS 

After all the analysis of the quasi-equatorial lensing 
in the Kerr spacetime, we are able to discuss the phe- 
nomenological relevance of the black hole spin. 

In Ref. 01 1 it was shown that all spherically sym- 
metric black holes produce the same patterns of rela- 
tivistic images. These patterns differ by the separations 
and the luminosity ratios between different relativistic 
images. The conclusion was that if strong field gravita- 
tional lensing will be caught by future VLBI experiments, 
it may provide a means to distinguish between different 
classes of black holes. 

In spherically symmetric black hole lensing, a generic 
source not aligned with the optical axis produces ex- 
tremely faint relativistic images. On the contrary, a point 
source perfectly aligned along the optical axis produces 
(theoretically) infinitely bright images. Actually the fi- 
nite source radius cuts off the real brightness of the im- 
ages. The relativistic images are maximally amplified 
altogether since all the caustic points of spherically sym- 
metric black holes lie on the optical axis. 

In Kerr lensing, the situation becomes radically differ- 
ent. The crucial fact is that the caustics no longer lie on 
the optical axis but drift throughout the trigonometric 
interval. Then if the source is close to one caustic point, 
it cannot be close to any other. The consequence is that 
only one image at a time can be enhanced, while all the 
others stay extremely faint. 

To clarify the situation, suppose we have a source 
aligned with the caustic point 73 ( | a | ) ■ Then the out- 
ermost relativistic image on the eastern side will be en- 
hanced. If we put the source on 7 3 (— |a|), then we only 
enhance the first relativistic image on the western side. If 
we put the source on 75 ( | a | ) , then the second relativistic 
image on the eastern side will be enhanced and so on. 

Rather than seeing an infinite series of relativistic im- 
ages on each side of the lens, we would observe only 
one enhanced relativistic image. It would be difficult to 
recognize a single image as a gravitational lensing phe- 
nomenon rather than any kind of environmental source 
around the black hole. Even if we managed, it would 
be quite tricky to extract information about the strong 
fields around the black hole from one single relativistic 
image. 

However, if the source is inside the caustic, two addi- 
tional images should appear, making easier to recognize 
their real nature of lensed images of the same source. 
Moreover, three images can be used to investigate the 
gravitational field around the black hole and put con- 
straints on its parameters. Alas, the additional images 
are missed in our quasi-equatorial approximation. In or- 
der to catch them, it is necessary to face the problem 
of Kerr lensing in its general form. A reliable treat- 
ment of non equatorial images, would complete our quasi- 
equatorial study and would make possible a detailed in- 
vestigation of strong field gravitational lensing for high 
values of the black hole spin. 
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By now, we have discussed the two extreme situations: 
spherically symmetric black holes (a = 0) and high spin 
black holes. It is interesting to estimate the value of the 
spin which separates the two regimes. We shall do it, 
referring to the black hole at the center of our Galaxy, 
assuming Dol — 8.5 kpc, Dls = 1 kpc. 

In order to consider a black hole as spherically symmet- 
ric, the caustic drift must be negligible when compared 
to the extension of the source. In fact, in this case, the 
source does not "see" different caustics but they behave 
roughly as they were all at the same point. For small a 
the 7^ scale linearly as —(k — l)ir — a(?7o + fc ) , with rjo 
and r]i numerical factors of order one. Then, the drift 
between two consecutive caustics will be negligible if 

|7 fc (o) - 7fc+2(a) mod 2vr| = 2 Vl a < (102) 

L>LS 

If we consider a source with radius Rs — lOi?©, we 
then find that a should be lower than 10 -10 . When a 
is greater than this value, the source will see only one 
caustic at a time. However, the caustic will be still seen 
as point-like, since the extension of the caustic scales as 
a 2 . Therefore, we will still have two enhanced images, 
corresponding to the two intersections of the fc-th caus- 
tic with the equatorial plane 7fc(a) and 7&(— a). This 
intermediate situation takes place as long as 

|7 fc (o) + (fc-l)7r| -|7fc(-o) + (fc-l)7r| = 

(£o+6fc)a 2 «^ (103) 

with £o an d £i of order one. 

For a lOi?© source, this requires a <C 10 -5 . Beyond 
this value, only one image at a time will be enhanced 
(together with an eventual additional pair of images if 
the source is inside the caustic), while all the others stay 
invisible. 

These estimates reveal that the phenomenology of 
spherically symmetric black holes is realistic only for 
black holes with tiny spin. Yet, as recalled in the in- 
troduction, the first estimates of the spin of the black 
hole at the center of our Galaxy push towards high val- 
ues [l4(- If these estimates are confirmed, then we are 
forced to include spin in any realistic treatment of strong 
field gravitational lensing for Sgr A*. 

X. SUMMARY 

In this paper we have explored the modifications to 
strong field limit gravitational lensing induced by the ro- 



tation of the central body, analyzing the quasi-equatorial 
null geodesies. 

The most apparent change is the formation of extended 
caustics which, for high angular momenta, can cover sev- 
eral Riemann folds. This situation is radically different 
from spherically symmetric black holes where the caustics 
are points aligned behind and in front of the lens. While 
for a = 0a source behind the lens is simultaneously close 
to all odd caustics and gives rise only to enhanced im- 
ages, for Kerr black holes the source can be close to one 
caustic at a time and thus produces only one enhanced 
image. 

As secondary interesting effects, we can also mention 
the asymmetry between images formed by photons wind- 
ing in the same sense of the black hole and photons wind- 
ing in the opposite sense, the latter appearing farther 
from the black hole. The magnification decreases with 
the spin, being higher for retrograde images. 

The study of quasi-equatorial Kerr gravitational lens- 
ing is very instructive and has allowed us to discover a 
great number of interesting features of spinning black 
holes. However, to address the phenomenology of the 
black hole at the center of our Galaxy and/or other black 
holes with deeper insight, further investigation is neces- 
sary. In fact, we need a punctual description of the caus- 
tic structure not limited to the equatorial plane. The 
existence of extended caustics suggests the formation of 
pairs of non-equatorial images, missed in our approxi- 
mation, which are of striking importance for the phe- 
nomenology. 

The quasi-equatorial lensing, studied in this work, then 
represents a first fundamental step to understand lens- 
ing by spinning black holes. However, the complexity of 
the problem requires a global approach in order to give 
correct and complete answers to all observational ques- 
tions. This settles as the main objective for future work 
on strong field gravitational lensing. 
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